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Abstract 

We present the results from a calculation of the full electroweak two-loop fermionic 
contributions to the effective leptonic mixing angle of the Z b oson, sin ^efF? 
the Standard Model. On-shell renormalization and analytic calculations are per- 
formed for the three-point vertex functions at zero external momenta, whereas ir- 
reducible three-point integrals for non-vanishing external momenta are evaluated 
semi-analytically applying two different methods. Comparisons with a previous cal- 
culation show complete agreement. 



1 Introduction 



Precision observablcs at the Z-boson resonance, measured at LEP 1 and SLC, together 
with the masses of the W boson and of the top quark, measured at LEP 2 and the Tevatron, 
constitute a set of high-energy quantites which in comparison with the Standard Model 
predictions aUow global precision analyses yielding indirect bounds on the Higgs-boson 
mass Mh [!]• In this context, one of the most important quantities with a very high 
sensitivity to Mh is the effective leptonic mixing angle, expressed as sin^ ^cflfj in the effective 
couplings of the Z to leptons (electrons, to be precise). Besides Mw, it is a key observable 
measured with high experimental accuracy, with expected further improvements in the 
future. The current experimental value is 0.23147 ± 0.00017 [1]; a linear electron-positron 
collider with GigaZ capabilities could even reach an accuracy of 1.3 x 10^^ [2, 3]. Therefore, 
very precise and reliable theoretical predictions for sin^ ^efr are required. 

The electroweak mixing angle in the effective leptonic vertex of the Z boson can be 
defined via the relation (see e.g. [4]) 

sin^^eff = Ul-Re^), (1) 
4 V 9aJ 

in terms of the dressed vector and axial vector couplings gv,A in the Z-lepton vertex. 
Starting from the Z mass Mz and the W mass Mw, sin^^eflf can be obtained in the 
following way, 

/ M^ \ 

sin^^eff = i^slr = 1^ \^^- J^J , k=1 + Ak, (2) 

with higher-order contributions collected in the quantity Ak. The W mass can be related 
to the precisely known Fermi constant Gp with the help of the W-Z interdependence, 

where Ar denotes the higher-order contributions. Recently, the full electroweak two-loop 
calculation for Ar has been completed, with both fermionic loops [5, 6] and purely bosonic 
loops [7] , and further improvement by three-loop contributions to the p-parameter [8] and 
to the S parameter [9] has been given. Moreover, the universal QCD corrections via the 
self energies [10] are known. 

For the k factor in (2), fewer higher-order terms have been evaluated as yet. Besides 
the universal contributions through the p, S parameters and the QCD corrections, a top 
mass expansion for the electroweak two- loop corrections was performed [11], and only 
recently the complete two-loop electroweak corrections of the fermionic type, i.e. with at 
least one closed fermion loop, were calculated [12]. In this paper we present an independent 
calculation of those fermionic two-loop contributions and perform a comparison with the 
result of [12]. 
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2 Electroweak two-loop contributions 



Expanding the dressed couplings in (1) according to gy^A = fJv^A y- + Qv^A + 9v^A + " " " j 
powers of a, one obtains the O (a^) contribution to sin^ ^eflf in the loop expansion 

sin^^eff = sin^ + sin^ ^2 + sin' ^2 + (4) 

to be 

(0) 

«in^« = -j^R^{s'^'-9T+£'{s^'-9'^'))- (8) 

Hence, besides the two-loop diagrams depicted schematically in Fig. 1, products of one- 
loop contributions have to be taken into account. The circles in the diagrams of Fig. 1 
denote the renormalized two- and three-point vertex functions, at the one-loop level in the 
reducible diagrams of Fig. la and lb, and at two-loop order in the irreducible diagrams of 
Fig. Ic and Id. 

In our computation we have used the on-shell renormalization scheme as described 
in [5]. In this scheme the real part of the diagram shown in Fig. Ic vanishes. The reducible 
diagrams of Fig. la and lb only contribute products of imaginary parts of one-loop func- 
tions, which can be easily computed. So the genuine two-loop task is the calculation of the 
irreducible Zii-vertex diagrams in Fig. Id. 

For the computation of the two-loop corrections with at least one closed fermion loop 
we are directed to the vertex diagrams shown in Fig. 2, one- loop diagrams with one- loop 
counter term insertions of the fermionic type, and the two-loop fermionic counter term. In 
the combination entering (5), the two-loop counter term for sin^ ^eflf in the on-shell scheme 
is given by 



-|- products of one loop counter terms. 

SZ'^^ and 5Z^^^^^ are the field renormalization constants for the photon-Z two-point 
function and for the left-/right-handed lepton field. The counter term for the on-shell 
quantity in (2) contains the two-loop counter terms for the Z- and VF-boson masses, 
(5M|^^^2)) ^nd products of one- loop counter terms. 
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Apart from the products of one-loop functions in (5) and the reducible diagram Fig. lb 
only the two-loop counter term (6) yields contributions with two closed fermion loops. 

The contributions with one closed fermion loop contain the generic two-loop diagrams 
depicted in Fig. 2. Their evaluation encounters a twofold task: the two-loop renormaliza- 
tion, and the computation of genuine two-loop vertex functions for non- vanishing external 
momenta. We have separated these two items by splitting the renormalized two-loop con- 




c) d) 
Figure 1: Classes of two loop diagrams 
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tribution to the Zii vertex function with the two external lepton momenta p±, 



Vl^'^\P') = 7. [9^ - 9'^!^ , - {P- +P^r , (8) 

into two UV-finite pieces according to 

= [rJ«(2)(o) + rf] + [rf^(^)(M|)-rf^(^)(o)]. (9) 

pZ«(2) ^p2-^ denotes the corresponding unrenormahzed Z££ vertex, and F^^ is the two-loop 
counter term, which is independent of P^. The first term on the right hand side of (9) 
contains the complete two-loop renormalization, but no genuine two-loop vertex diagrams 
since in absence of external momenta they reduce to simpler vacuum integrals. All the 
genuine two-loop vertex diagrams appear as subtracted quantities in the second term on 
the right hand side of (9). 




c) d) 
Figure 2: Irreducible two-loop vertex diagrams containing fermion loops. 
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3 Outline of the calculation 

In order to calculate the irreducible vertex diagrams of Fig. 2, various techniques are 
applied. Diagrams containing a one-loop self-energy insertion are calculated with the help 
of dispersion relations. This results in one-dimensional integrals, which are computed 
numerically. 

The diagrams with fermion triangles, Fig. 2d, require a careful treatment of 75 in 
D dimensions. In [5] it was shown that naive dimensional regularization with a four- 
dimensional treatment of the e-tensor from traces involving 75 can be applied for this 
special case. Putting such traces to zero yields a UV-finite difference, which vanishes for a 
fermion generation with equal masses. Since the light fermion masses have been neglected 
in our approach, a non-vanishing contribution occurs only from the third generation. In 
the numerical discussion of section 4 this difference is labeled as the "75" contribution. 

Differently from the charged-current case treated in [5] , the neutral-current vertex con- 
sidered here involves another subtlety originating from photons in the internal lines of 
Fig. 2d, which give rise to collinear divergences in the case of massless electrons. This 
would produce unphysical finite extra shifts in combination with the {D — 4) terms from 
the 75 traces. [For massive internal bosons those terms disappear in the limit D — > 4 
because the second loop integration does not introduce another UV divergence.] In order 
to prevent this kind of problems, we keep the physical mass of the external electrons in 
the dangerous diagrams. Moreover, also the masses of the internal b quarks and r leptons 
have been kept. The resulting two-loop vertex functions are further evaluated using the 
methods described in [13]. In the end we are left with up to four-dimensional integrals, 
which are calculated numerically. 



In order to get the renormalized vertex (respectively, the contribution to sin^ ^eff) from the 
two-loop vertex at zero momentum, the two-loop rcnormalization constants in (6) and (7) 
are needed. In the on-shell scheme [5] they arc given by 



5M^(,) = Re(E;^)(M^))-5Z(T)5M^(i)+Im(s;r)\M^))lm(s;[)(M^^^ (10) 



3.1 Computation of ff^^(2)(o) 



5Ml 



Z{2) 



Re (Egf 



(M|))-5Z(if5M|(,)+Im 



(Egf(M|))lm(Egf(M|)) 




(11) 



(12) 
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involving the bosonic self energies and field renormalization constants, and 

5Z'^^ = -Sj,«(0), (14) 

from the two-loop leptonic self energies according to the decomposition 

H'ip) =^u;iE^^(p2)+ ^u;flE^^(/) . (15) 

E' ip^) denotes the derivative of a self-energy E (p^) with respect to p^. 

The computation of the renormalized vertex at vanishing external momentum is done in 
the following steps. The Feynman diagrams are generated with the package FeynArts [14]. 
The program TwoCalc [15] is applied to reduce the amplitudes to standard integrals. 
The resulting scalar one-loop integrals and two-loop vacuum integrals are calculated using 
analytic results [16, 17]. The two-loop self-energy functions with non- vanishing external 
momenta appearing in the renormalization constants (10)-(12) are obtained with the help 
of one-dimensional integral representations [18]. 

3.2 Diagrams with self-energy insertions 

The strategy to compute vertex diagrams with self-energy insertions for the subtracted 
vertex r^^^(Mf ) — r^^^(O) will be illustrated in terms of Fig. 2b with a diagonal Z self- 
energy as an example (the other diagrams are calculated in an analogous way). The 
corresponding vertex diagram has the structure 

/I 1 ( 1 

d^g Yigv-9al5) ^_ ^_ li^igv-Qal^) l" (^^-^a75) ^^2^71^2)2 ' (1^) 

with = (j9_ +P+)^- The renormalized one-loop self-energy can be decomposed into its 
transverse and longitudinal part, 

Sp.(?) = St(0 (^P. - ^) + S.(0 (^) • (17) 

It can be easily verified that for on-shell leptons the part proportional to gpg^r yields a term 
independent of P^, hence it drops out in the difference r^^^(M|) — r^^^(O), and only the 
transverse part has to be taken into account. In the on-shell-scheme it reads 

ET(g') = ST(g')-ReST(M^)-(g2-M|)ReE^(M|) (18) 
= Eriq') - T,t{MI) - {q' - M|) E^(M|) 
+i ImEr(M|) + i {q^ - M|) ImE^(M|) . 

Inserting the imaginary parts into (16), one gets a constant times one-loop functions. For 
the rest, we use the dispersion relation 

, I f°° , ImST(s) , , 

Et = - / ds f^, 19 

TT Jo s-q^ -le 
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yielding 



{q^-Mlf ImET(s) 1 



Inserting this into (16), the factor (g^ — M|)^ cancels the corresponding term in the de- 
nominator, and we end up with a one-dimensional integral over a one-loop three-point 
function with one variable mass s, multiplied by the weight-factor ImS (s) /{s — M|)^. 
The subtraction at = ensures that the expressions are UV finite. Finally we are left 
with integrals of the type 



His , = / 

Jo (s - M2 - ief Jo 



f is) r / (^) - / (M^) - (. - M^)r (M^) 



s - M2 - iey Jo {s-M^- ie) 

^ j'""" f{M^) + {s-M^)f{M^) 
{s-M^- ief 







The second term in (21) can be computed analytically. The other terms in (21) arc free 
of singularities and can be integrated numerically. We have used the CUBA library [19] 
which allows a very fast and precise evaluation of the integrals. 



3.3 Diagrams with fermion triangles 

For the computation of the diagrams containing fermion triangles (Fig. 2d) we have basi- 
cally adopted the numerical methods described in [13]. The mass of the external electron 
lines is set to zero except for the cases with internal and 77, where rrie is used to 
regularize infrared and coUinear divergences. In the cases Z7 and 77 we have also kept 
all the internal fermion masses (6, t,r), of the third generation, while in the cases WW 
and Z Z all fcrmions are treated as massless, except the top quark. 

The resulting expression for the contribution to sin^ 6^ in (1) can be written as a sum 
of terms from tensor integrals of the family V'^^^ (the families of two-loop self-energies 
and vertices are collected in the appendix). In order to simplify as much as possible the 
tensorial structure, we perform a simple reduction of the type 

2q-p _ 1 1 p2 - 

(g2 + rn2) [(g + p)2 + M2] ~ g2 + m2 ~ (q + pf + M2 ~ (g2 + rn?) [(^ + p)2 + m^] ' 

which does not introduce any new denominator and therefore any spurious singularity. At 
this point we have a sum of tensor integrals up to rank 3 belonging to four vertex families 
^yi2i^ yi3i^ y22i^ ^.231)^ together with three self-energy types {S^^^, S^'^^, S^^^), and some 
one-loop diagrams {B and C functions). 
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For all these tensor integrals, Feynman parameters are introduced. The parameter in- 
tegrals after integration over the loop momenta are manipulated in order to obtain smooth 
integrands to be computed numerically. The methods used are those described in [13] to- 
gether with some improvement aimed to increase the numerical stability in specific regions 
of the phase-space (some new methods have been also developped for internal checks). In 
this context, it is worth mentioning the discovery of a new kind of algorithm for the com- 
putation of one-loop diagrams, which has been presented in [20]. The basic idea is that, 
given a polynomial V in the Feynman parameters x — {xi,X2, ...)*, a constant B and a 
column vector V with 

V{x)^Q{x)+B, V' Q{x) = -Q{x) , (22) 

it can be easily proved that the following relation holds {(3 > 0), 



J dxV{xY = J dx dyy^-^ [Q{x)y + B]'' . 



(23) 



Since any one-loop diagram can be expressed as an integral of a quadratic polynomial to 
a negative power, the previous formula can be applied with /i < 0. After an integration 
by parts, which cancels d^-, we can perform the integration in y. Then the procedure can 
be iterated with appropriate values of j3 to obtain smooth integrals. For example, for the 
scalar one-loop three-point function we obtain {e — A: — D) 



where 



V{x) = x^ Hx + 2K^x + L, B = L~K^H~^K, X = -H"^ K, 

V[0]{x,) = V{l,xi), V[l]{x,) = Vix^,x,), V[2]{x,) = V{x,,0), 

ao = l-Xi, ai = Xi-X2, a2 = X2, X = (Xi,X2)*. (25) 

The matrix H, the cohimn K, and the constant L are related to the physical quantities in 
the vertex diagram [Fig. 3] as follows (j, I — 1, 2), 

Hij = -pi-Pj, Kj = ]^{k]-k]_^ + m]^^-m]), L = ml - ie, 

ko = 0, ki^pi, k2^pi+P2- (26) 

This representation for Co (which can be easily generahsed for tensor integrals) is used 
to compute those two-loop diagrams that can be expressed as an integral of a one-loop 
three-point function, as in V'^'^^ and V"^^^ (see [13] for explicit expressions). In such cases 
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p 




Pi 



Figure 3: One-loop three-point function. The momenta are understood as incoming. 



the masses and momenta of the C functions depend on Feynman parameters, which have to 
be integrated over. It is therefore crucial that the zeros of the denominator — B) 

in (24) are compensated by the logarithm. 

Another important feature of the calculation is the infrared problem. Although the 
diagrams of Fig. 2d are IR finite, some individual integrals arising from the reduction in 
the Z 7 case are divergent. A typical case is the scalar part of the following diagram, which 
has the structure 



For all the IR divergent integrals in the decomposition we have extracted the infrared one- 
loop functions and verified analytically that the sum of them cancels. The remaining finite 
terms are then again expressed in terms of smooth integrals and computed numerically. 
Details of this treatment of IR divergences will be given elsewhere [21]. 

The algebraic handling and the numerical evaluation has been done in parallel by two 
independent computations. For the numerical integration the NAG library DOIGDF [22] 
was used in one case and the CUBA library [19] in the other one. According to the 
features of the integrator, the contribution to the final result are summed up either in 
a single integral or split into a sum of several terms. In the second case the total error 
has been computed as the sum of the errors of the different integrals. Thanks to the use 
of more than one method for each diagram type, it has been possible to achieve a good 
numerical precision and make a severe cross-check of the calculations. 




finite terms. (27) 
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4 Results 



The input parameters for the evaluation of the final result are put together in Tab. 1. They 
are chosen in accordance with [12] in order to make an immediate comparison possible. Mw 
and Mz are the experimental values of the W- and Z-boson masses [23], which are the on- 
shell masses. They have to be converted to the values in the pole mass scheme [5], labeled 
as Mw and Mz- which are used internally for the calculation. These quantities are related 
via Mw,z = Mw,z + r^^/(2 Mw,z)- For Vz the experimental value (Tab. 1) and for Vw 
the theoretical value has been used, i.e. Tw = 3 G^My^/ (2\/27r) (1 + 2as (M^) / (37r)) 
with sufficient accuracy. 



parameter value 

~M^ 80.426 GcV 

Mz 91.1876 GeV 

Vz 2.4952 GeV 

mt 178.0 GeV 

Aa(M|) 0.05907 

as{Ml) 0.117 

1.16637 x 10-^ 

'Mw 80.3986 GeV 

Mz 91.1535 GeV 



Table 1: Input parameters entering our computation. Mw and Mz are the experimental values 
of the W- and Z-boson masses, whereas Mw and M z are the calculated quantities in the pole 
mass scheme. 

The results are given for A/t, eq. (2), and are listed in Tab. 2 and Tab. 3 for various 
masses of the Higgs boson. Tab. 2 contains also the one-loop result for comparison and the 
corresponding results obtained in [12]. A finite 6-quark mass has been kept in the one- loop 
result. The errors on our two- loop-results (in brackets) are due to the uncertainty from 
the numerical integration. Our results are in full agreement with those given in [12] (last 
column of Tab. 2). 

In Tab. 3 the various parts of the two-loop result are shown. Large cancellations between 
the part containing two fermion loops and the part containing only one fermion loop occur. 
Moreover, one can see that, depending on the value of Mh, the first term of (9) (column 4) 
is about 15-20 times larger than the second one (column 5). Hence, the complete result for 
the renormalized Zll vertex at = M| can be well approximated by the much simpler 
expression at = 0. The "75" part is also very well approximated by its value at = 0, 
which is 0.280 x 10""^. The "75" part does not contain any terms proportional to or 
m^; for illustration it is plotted in Fig. 4, at P^ = as a function of the top mass. 
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Mh [GeV] O {a) x 10"^ 


O (a^) x 10-^ 


0(q;2) x 10- 


4ri oi 
[12J 


100 


438.937 


-0.637(1) 


-0.63 




200 


419.599 


-2.165(1) 


-2.16 




600 


379.560 


-5.012(1) 


-5.01 




1000 


358.619 


-4.737(1) 


-4.73 




Table 2: 


Two-loop result for A/c in comparison with the one-loop result and the result in 


[12]. 


Mh 
[GeV] 


2 ferm. loops red. 

xlO~^ xlO-^ 


f (0) r (M|) - r (0) 75 
xlO"^ xlO"^ xlO"^ 




100 


13.758 -0.722 


-14.903 0.959(1) 


0.271 




200 


13.758 -0.688 


-16.465 0.959(1) 


0.271 




600 


13.758 -0.501 


-19.499 0.959(1) 


0.271 




1000 


13.758 -0.386 


-19.339 0.959(1) 


0.271 





Table 3: Breakdown of the two-loop result. The column "2 ferm. loops" gives all contributions 
with two closed fermion loops, whereas the other columns only contain results with one closed 
fermion loop, "red." refers to the reducible contribution of the diagrams depicted in Fig. la 
and Fig. lb plus the product of one-loop contributions in (5). F (0) is the first term of (9) and 
F (M|) — F (0) the second one. In the last column the "75 "-contributions are given (see text). 




200 400 600 800 1000 

mt [GeV] 



Figure 4: "75" contribution at = as a functions of mj. 



In conclusion, we have evaluated the electroweak 2-loop corrections to sin^ ^eff with at least 
one closed fcrmion loop. Methods to calculate the appearing two-loop vertex functions have 
been described. A discussion of the various individual parts of the two-loop result was given 
and agreement with [12] was found. 
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12 



A Two-loop self-energy and vertex families 
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